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Abstract 

These notes give a few examples illustrating how continuum models can be derived from special 
limits of discrete models. Only the simplest cases are considered, illustrating some of the most 
basic ideas. These techniques are useful because continuum models are often much easier to deal 
with than discrete models with very many variables, both conceptually and computationally. 
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1 Introduction. 

Continuum approximations are useful in describing discrete systems with a large number of degrees 
of freedom. In general, a continuum approximation will not describe all possible solutions of the 
discrete system, but some special class that will depend on the approximations and assumptions 
made in deriving the continuum model. Whether or not the approximation is useful in describing 
a particular situation, will depend on the appropriate approximations being made. The most 
successful models arise in situations where most solutions of the discrete model evolve rapidly in 
time towards configurations where the assumptions behind the continuum model apply. 

The basic step in obtaining a continuum model from a discrete system, is to identify some basic 
configuration (solution of the discrete model) that can be described by a few parameters. Then one 
assumes that the full solution of the system can be described, near every point in space and at every 
time, by this configuration — for some value of the parameters. The parameters are then assumed 
to vary in space and time, but on scales (macro-scales) that are much larger than the ones associated 
with the basic configuration (micro-scales). Then one attempts to derive equations describing the 
evolution of these parameters in the macro-scales, thus averaging out of the problem the micro- 
scales. There is a close connection between this approach, and the"quasi-equilibrium" approximations 
that are often invoked to "close" continuum sets of equations derived using conservation laws. 

For example, when deriving the equations for Gas Dynamics in Statistical Mechanics, it is assumed 
that the local particle interactions rapidly exchange energy and momentum between the molecules 
— so that the local probability distributions for velocities take a standard form (equivalent to local 
thermodynamic equilibrium). What exactly makes these assumptions work (in terms of properties 
of the governing, micro-scale, equations) is rather poorly understood. But that they work rather 
well cannot be denied. In these notes we will consider examples that are rather simpler than these 
ones, however, where the "local configurations" tend to be rather trivial. 

2 Wave Equations from Mass-Spring Systems. 

Longitudinal Motion. 

Consider an array of bodies/particles, connected by springs, and restricted^ to move on a straight 
line. Let the positions of the bodies be given by Xn = Xn{t), with n = 0, ±1, ±2, . . ., and let Mn 
be the mass of the n*'* particle. Furthermore, let the force law for the spring between particles 
n and n + 1 be given by: force = /„_,_i(Aa:), where Ax is the distance between the particles, and 
is positive when the spring is under tension.^ 

If there are no other forces involved (e.g. no friction), the governing equations for the system are: 



for n = 0, ±1, ±2, . . . The simplest solution for this system of equations is equilibrium. In this 
case all the accelerations vanish, so that the particle positions are given by the series of algebraic 

^By some device: say the bodies are sliding inside a hollow tube. 



^If the spring obeys Hooke's law, then /„_,_i(Aa;) = \^Ax — L„_|_i j , where > and i„-|-i > are the 

spring constant and equilibrium length, respectively. 




(2.1) 
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equations 

= /„+ 1 {Xn+l - Xn) - fn-l{Xn " Xn-l) ■ (2.2) 

This is the basic configuration (solution) that we will use in obtaining a continuum approximation. 
Note that this is a one paranneter fannily: if the forces are monotone functions of the displacements 
Arr, then once any one of them is given, the others follow from (2.2). 

Before proceeding any further, it is a good idea to non-dimensionalize the equations. We will 
assume that: 

A. All the springs are roughly similar, so that we can talk of a typical spring force /, and a typical 
spring length L. Thus we can write 

fn^iM = f Fn+^, {^) , (2.3) 

where is a non-dimensional mathematical function, of 0(1) size, and with 0(1) deriva- 

tives. A further assumption is that F^_^_i changes slowly with n, so that two nearby springs 
are nearly equal. Mathematically, this is specified by stating that: 

F„+i(ry)=F(e(n + l/2),ry), (2.4) 
where < e ^ 1, and F is a "nice" (mathematical) function of its two variables. 

B. All the particles have roughly the same mass m, and their masses change slowly 
with n, so that we can write: 

Mn = mM{en), (2.5) 
where M is a nice mathematical function, with 0(1) size, and with 0(1) derivatives. 

Remark 2.1 Why do we need these assumptions? This has to do with the questions of validity, 
discussed in the introduction. Suppose that these hypothesis are violated, with the masses and 
springs jumping wildly in characteristics. Then the basic configuration described by (2.2) will still 
be a solution. However, as soon as there is any significant motion, neighboring parts of the chain 
will respond very differently, and the solution will move away from the local equilibrium implied by 
(2.2). There is no known method to, generically, deal with these sort of problems — which turn out 
to be very important: see remark 2.2. 

From the assumptions in A and B above, we see that: 



Changes in the nnass-spring systenn occur over length scales 



I = L/e. (2.6) 



Using this scale to non-dimensionalize space, namely: | Xn = iXn, | and a yet to be specified time 
scale T to non-dimensionalize time, namely: t = tT,] the equations become: 

A and B above also imply that, for the solution in (2.2), the inter-particle distance Xn+i — Xn varies 
slowly — an 0(e) fractional amount per step in n. Thus we propose solutions for (2.7) of the form: 

Xn{t) = X(s„, t) , where Sn = ne, (2.8) 

and X = X{s,t) is some smooth function of its arguments. 
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Substituting (2.8) into (2.7), and using (2.4) and (2.5), we obtain 




(2.9) 



Here we have used that: 




= -^X{s + \e,t) + 0{e^) and 



d_ 



■X{s 



\e,t) + 0{^), 



with a similar formula applying to the difference ~ ^n-^- 

Equation (2.9) suggests that we should take 



m L 



(2.10) 



T = 



for the un-specified time scale in (2.7). Then equation (2.9) leads to the continuum limit ap- 
proximations (valid for < e <C 1) 



The mass-spring system introduced in equation (2.1) can be thought of as a simple model for an 
elastic rod under (only) longitudinal forces. Then we see that (2.11) is a model (nonlinear wave) 
equation for the longitudinal vibrations of an elastic rod, with ,s a lagrangian coordinate 
for the points in the rod, M = M{s) the mass density along the rod, and X giving the position of 
the point s as a function of time, and F a function characterizing the elastic response of the rod. 
Of course, in practice F must be obtained from laboratory measurements. 

Remark 2.2 The way in which the equations for nonlinear elasticity can be derived for a crystalline 
solid is not too different? from the derivation of the wave equation (2.11) for longitudinal vibrations. 
Then a very important question arises (see first paragraph in section 1): What important behaviors 
are missed due to the assumptions in the derivation? How can they be modeled? In particular, 
what happens if there are "defects" in the crystal structure (see remark 2.1)? These are all very 
important, and open, problems of current research interest. 

Example 2.1 Uniform Rod. 

// all the springs and all the particles are equal, then we can take M = 1 and F is independent of 
s. Furthermore, if we take L to be the (common) equilibrium length of the springs, we then have 



where c^ = c^^q) = dF/dr]{iq) > 0, and F{1) = (equilibrium length). The unperturbed "rod" cor- 
responds to X = s, while X = as corresponds to the rod under uniform tension (a> 1), or com- 
pression (a<l). Also, note that c is a (non-dimensional) speed — the speed at which elastic 
disturbances along the rod propagate: i.e. the sound speed. 




(2.11) 




(2.12) 



^At least qualitatively, though it is technically far more challenging. 
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Example 2.2 Small Disturbances. 

Consider a uniform rod in a situation 'where tfie departures from uniform equilibrium are small. 
That is dX/ds ~ a, where a is a constant. Then equation (2.12) can be approximated by the 
linear wave equation 

Xtt = c^X,,, (2.13) 
where c = c{a) is a constant. The general solution to this equation has the form 

X = g{s-cT) + h{s + cT), (2.14) 

where g and h are arbitrary functions. This solution clearly shows that c is the wave propagation 
velocity. 

Remark 2.3 Fast vibrations. 

The vibration frequency for a typical mass m, attached to a typical spring in the chain, is: 



^ = = - ■ (2-15) 

This corresponds to a time scale much shorter than the one involved in the solution in (2.8-2.11). 
What role do the nnotions in these scales play in the behavior of the solutions of (2.1), under the 
assumptions made earlier in A and B ? 

For real crystal lattices, which are definitely not one dimensional (as the one in (2.1)) these fast time 
scales correspond to thermal energy (energy stored in the local vibrations of the atoms, relative to 
their equilibrium positions). It is believed that the nonlinearities in the lattice act so as to randomize 
these vibrations, so that the energy they contain propagates as heat (diffuses). In one dimension, 
however, this does not generally happen, with the vibrations remaining coherent enough to propagate 
with a strong wave component. The actual processes involved are very poorly understood, and the 
statements just made result, mainly, from numerical experiments with nonlinear lattices. 



Just to be a bit more precise: consider the situation where all the masses are equal 
for all n, and all the springs are equal and satisfy Hooke's law (linear elasticity): 

(Ax 



Mn = m 



f^^{Ax) = k{Ax -L) = f^—-lj, (2.16) 

where k is the spring constant, L is the equilibrium length, and f = k L.^ Then equation (2.1) takes 
the form 

— .T„ = io'^ {Xn+1 - 2Xn + Xn-l) , (2.17) 

where uj is as in (2.15). Because this system is linear, we can write its general solution as a linear 
superposition of eigenmodes, which are solutions of the forrn^ 

Xn = exp{iKn — iat) , where a = ±2 to sin ^— ^ and — oo<k<oo is a constant. (2.18) 

These must he added to an eqii/idiJjrmm solution Xn = a L n = Sn, where a > is a constant. 
"Check that these are solutions. 
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Relative to the mean position Sn along the lattice, each solution in (2.18) can he written as 

Xn = exp(i Sn -iat) . 

Q.L 

Thus we see that it represents a wave of wavelength X = 27raL/K, and speed 

aLa 2aLco . f k\ 2c . f k\ 

propagating along the lattice — where c = aLuj is a speed. Note that the speed of propagation is a 
function of the wave-length — this phenomenon is know by the name o/ dispersion. We also note 
that the maximum frequency these eigenmodes can have is a = 2io, and corresponds to wavelengths 
of the order of the lattice separation.^ 

In the case of equations (2.16 - 2.17) there is no intrinsic e in the equations: it must arise from the 
initial conditions. That is to say: assume that the wavelength I with which the lattice is excited is 
much larger than the lattice equilibrium separation L, i.e. £^ L, with e = L/l. This corresponds to 
solutions (2.18) with k small. In this long wave limit we see that (2.19) implies that the solutions 
have the same wave speed c^ = ±c. This corresponds to the situation in (2.13 - 2.14). 

It is clear that, in the linear lattice situation described above, we cannot dismiss the fast vibration 
excitations (with frequencies of the order of co) as constituting some sort of energy "bath" to be 
interpreted as heat. The energy in these vibrations propagates as waves through the media, with 
speeds which are of the same order of magnitude as the sound waves equation (2.13) describes. 
Before the advent of computers it was believed that nonlinearity would destroy the coherence of 
these fast vibrations. Numerical experiments, however, have shown that this is not (generally) true 
for one dimensional lattices,^ though it seems to be true in higher dimensions. Exactly why, and 
how, this happens is a subject of some current interest. 

Transversal Motion. 

We consider now a slightly different situation, in which the masses are allowed to move only in 
the direction perpendicular to the x axis. To be precise: consider a sequence of nnasses M„ in 
the plane, whose x coordinates are given by Xn = nL. Each mass is restricted to move only in the 
orthogonal coordinate direction, with y„ = giving its y position. The masses are connected by 

springs, with (Ar„_,_i ) the force law, where Ar„_,_i = \J L'^ + [Un+i — UnY distance between 

masses. Assuming that there are no other forces involved, the governing equations for the system 
are: ^ 

d Un+l ~ Un J. If. N, Vn ~ Un-l r / \ \ /r, on\ 

"^"^^^ = - ' (2.20) 



"+2 ""2 



for n = 0, ±1, ±2, . . . (you should convince yourself that this is the case). 

The simplest solution for this system of equations is equilibrium, with all the masses lined up 
horizontally | yn+i = Vn, \ so that all the accelerations vanish. Again, one can use this (one 
parameter) family of solutions to obtain a continuum approximation for the system in (2.20) — 
under the same assumptions earlier in A and B. 



^The reason for the 2 relative to (2.15) is that the masses are coupled, and not attached to a single spring. 

^The first observation of this general phenomena was reported by E. Fermi, J. Pasta and S. Ulam, in 1955: Studies 
of Non Linear Problems, Los Alamos Report LA-1940 (1955), pp. 978-988 in Collected Papers of Enrico Fermi. II, 
The University of Chicago Press, Chicago, (1965). 
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Remark 2.4 Stability of the Equilibrium Solutions. 

It should be intuitively obvious that the equilibrium solutions described above will be stable only if the 
equilibriunn lengths of the springs C^j^i are smaller than the horizontal separation L between the masses, 
namely: < L. This so that none of the springs is under compression in the solution, since any 

mass in a situation where its springs are under compression will easily "pop " out of alignment with 
the others — see example 2.3. 

Introduce now the non-dimensional variables Y = ey/L, X = ex/L (note that, since Xn = nL, in 
fact X plays here the same role that s played in the prior derivation^), and T = t/r, where r is as 
in (2.10). Then the continuum limit for the equations in (2.20) is given by 

a^r d ( F{x, s) dY\ 
'^^^^^df^ = dx[-^—dx) ^^-^^^ 

where Y = Y{X,T) and 

S 



( byV 

1 + 



I The derivation of this equation is left as an exercise to the reader. 



The mass-spring system introduced in (2.20) can be thought of as a simple model for an elastic 
string restricted to move in the transversal direction only. Then we see that (2.21) is a model 
(nonlinear wave) equation for the transversal vibrations of a string, where X is the 
longitudinal coordinate along the string position, Y is the transversal coordinate, M = M(X) is 
the mass density along the string, and F = F{X, S) describes the elastic properties of the string.^ In 
the non-dimensional coordinates, the (local) equilibrium length for the string is given hj ei = C/L. 
That is, the elastic forces vanish for this length: 

F{X,ei{X)) = , where < 1 (for stability, see remark 2.4). (2.22) 



We also assume that 



Example 2.3 Uniform String with Small Disturbances. 

Consider now a uniform string (neither M , nor F, depend on X) in a situation where the departures 
from equilibrium are small (dY/dX is small). 

For a uniform string we can assume M = 1, and F is independent of X. Thus equation (2.21) 
reduces to 

d^^^fF{SldY\ 

dT^ dX { S dx) ' ^ ■ ^ 

Next, for small disturbances we have <S ~ 1, and (2.23) can be approximated by the linear wave 
equation 

Ytt = c:'Yxx, (2.24) 
wiiere = F{1) is a constant (see equations (2.13 - 2.14). 



^The coordinate s is simply a label for the masses. Since in this case the masses do not move horizontally, X can 
be used as the label. 

^Notice that <S is the local stretching of the string, due to its inclination relative to the horizontal position (actual 
length divided by horizontal length). 
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Notice how the stability condition < 1 in (2.22) guarantees that c^ > in (2.23). If this were not 
the case, instead of the linear wave equation, the linearized equation would have been of the form 

Frr + c?'Fxx = 0, (2.25) 

with d > 0. This is Laplace Equation, which is ill-posed as an evolution in time problem. 

To see this, it is enough to notice that (2.25) has the following solutions: 

Y = e'^l'^l* mn{kX) , for any - oo < k < oo . (2.26) 

These solutions grow arbitrarily fast in time, the fastest the shortest the wave-length (\k\ larger). 
This is just the mathematical form of the obvious physical fact that a straight string (with no bending 
strength) is not a very stable object when under compression. 



General Motion: Strings and Rods. 

If no restrictions to longitudinal (as in (2.1)) or transversal (as in (2.20)) motion are imposed on the 
mass-spring chain, then (in the continuum limit) general equations including both longitudinal and 
transversal modes of vibration for a string are obtained. Since strings have no bending strength, 
these equations will be well behaved only as long as the string is under tension everywhere. 

Bending strength is easily incorporated into the mass-spring chain model. Basically, what we need 
to do is to incorporate, at the location of each mass point, a bending spring. These springs apply 
a torque when their ends are bent, and will exert a force when-ever the chain is not straight. The 
continuum limit of a model like this will be equations describing the vibrations of a rod. 

We will not develop these model equations here. 

3 Torsion Coupled Pendulums: Sine-Gordon Equation. 

Consider an horizontal axle A, of total length i, suspended at its ends by "friction less" bearings. Along 
this axle, at equally spaced intervals, there are N equal pendulums. Each pendulum consists of a rigid 
rod, attached perpendicularly to the axle, with a mass at the end. When at rest, all the pendulums 
point down the vertical. We now make the following assumptions and approximations: 

M 

• 1. Each pendulum has a mass — . The distance from its center of mass to the axle center is L. 

• 2. The axle A is free to rotate, and we can ignore any frictional forces (i.e.: they are small). In 

fact, the only forces that we will consider are gravity, and the torsional forces induced on the 
axle when the pendulums are not all aligned. 

• 3. Any deformations to the axle and rod shapes are small enough that we can ignore them. Thus 

the axle and rod are assumed straight at all times. 

• 4. The mass of the axle is small compared to M , so we ignore it (this assumption is not strictly 

needed, but we make it to keep matters simple). 

Our aim is to produce a continuum approximation for this system, as AT ^ oo, with everything else fixed. 



Discrete to Continuum Modeling. MIT, March, 2001 — Rosales. 



9 



Each one of the pendulums can be characterized by the angle 6n = On{t) that its suspending 
rod makes with the vertical direction. Each pendulum is then subject to three forces: 

(a) Gravity, for which only the component perpendicular to the pendulum rod is considered.^ 

(b) Axle torsional force due to the twist 6n+i — On- This couples each pendulum to the next one. 

(c) Axle torsional force due to the twist 9n — On-i- This couples each pendulum to the prior one. 

We will assume that the amount of twist per unit length in the axle is small, so that Hooke's law applies. 
Remark 3.1 Hooke's Law for Torsional Forces. 

In the Hooke's law regime, for a given fixed bar, the torque generated is directly proportional to the 
angle of twist, and inversely proportional to the distance over which the twist occurs. 
To be specific: in the problem here, imagine that a section of length Ai of the axle has been twisted 
by an amount (angle) ^. Then, ifT is the torque generated by this twist, one can write 



K 



M 



(3.1) 



where k is a constant that depends on the axle material and the area of its cross-section — assume 
that the axle is an homogeneous cylinder. The dimensions of k are given by: 

. , mass X lenqth^ force x area 

M = — — 9 ^ = ■ 3.2 

time X angle angle 

This torque then translates onto a tangential force of magnitude F = T/L, on a mass attached to 
the axle at a distance L. The sign of the force is such that it opposes the twist. 

Let us now go back to our problem, and write the equations of motion for the pendulums. We 
will assume that: 

i 

• The horizontal separation between pendulums is . 

^ ^ N + 1 

£ 

• The first and last pendulum are at a distance —— from the respective ends of the axle. 

2(A^ + 1) ^ 

The tangential force (perpendicular to the pendulum rod) due to gravity on each of the masses is 

Fg = -^MgsinOn, where n = l,...,N. (3.3) 

For any two successive masses, there is also a torque whenever 9n 7^ On+i- This is generated by the 
twist in the axle, of magnitude 9n+i — 9n, over the segment of length £/ (N + 1) connecting the two 
rods. Thus each of the masses experiences a force (equal in magnitude and opposite in sign) 

Ft = ±{N + 1)^{0^+,-0^), (3.4) 

where the signs are such that the forces tend to make 9n = 9n+i. Putting all this together, we obtain 
the following set of equations for the angles: 

h""-^ = -i^S-^. + ^^^C-".), (3.5) 



^The component along the rod is balanced by the rod itself, which we approximate as being rigid. 
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LmL^ = -1- Mg smOn 
N dt^ N ^ 

+ - iPn+l - On) - ^ {On " ^n-l) , (3.6) 

for n = 2, . . . , AT — 1, and 

iMii^ = -^Mssme«-(^^^(e„ (3.7) 
These are the equations for N torsion coupled equal pendulums. 

Remark 3.2 To check that the signs for the torsion forces selected in these equations are correct, 
take the difference between the n^^ and (n + 1)^^ equation. Then you should see that the torsion 
force (due to the portion of the axle connecting the n^^^ and (n + 1)^^ pendulums) is acting so as to 
make the angles equal. 

Remark 3.3 Note that the equations for the first and last angle are different, because the first and 
last pendulum experience a torsion force from only one side. How would you modify these 
equations to account for having one (or both) ends of the axle fixed? 

Continuum Limit. 

Now we consider the continuum limit, in which we let AT ^ oo and assume that the angle 
can be written in the form: 

en{t) = e{xn,t), (3.8) 

n + - 

where 9 = d{x,t) is a "nice" function (with derivatives) and Xn = ^ ^ ^ is the position of the 
pendulum along the axle. In particular, note that: 

£ 

Ax = Xn+l -Xn = Y ■ (3-9) 

Take equation (3.6), and multiply it by N/i. Then we obtain 

r d'^^n . . N(N + 1)k X 
dt^ = -pgsm9n + -pj^ [6n+l - ^On + On-l) , 



where 



p = M/£ 



is the mass density per unit length in the N ^ oo limit. Using equation 



(3.9), this can be written in the form: 

pL^ = -P9^^-0n + ^ . (3.10) 

From equation (3.8) we see that — in the limit N ^ oo (where A ^ 0) — we have: 



(Ax)2 dx^ 
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Thus, finally, we obtain (for the continuum limit) the nonlinear wave equation (the "Sine— Gordon" 
equation): 

9u-c^9xx = -cousin 9, (3.11) 



where 



is the pendulum angular frequency, and 



K 



is a wave propagation speed 



(check that the dimensions are correct). 
Remark 3.4 Boundary Conditions. 

What happens with the first (3.5) and last (3.7) equations in the limit N oo? 

As above, multiply (3.5) by l/l. Then the equation becomes: 

pLd^9i_ pg {N + 1)k pg k 92 - 9i 

AT ^ - " AT + - - ^' + JE -K^ 

Thus, as N ^ oo one obtains 



9x{0J) = 0. 



This is just the statement that there are no torsion forces at the x = end (since the axle is free to 
rotate there). Similarly, one obtains: 



9x{Lt) = 0, 



at the other end of the axle. How would these boundary conditions be modified if the axle 
where fixed at one (or both) ends? 



Kinks and Breathers for the Sine Gordon Equation. 

Equation (3.11), whose non-dimensional form is 



9u — 9rr 



sin ^ . 



(3.12) 



has a rather interesting history. Its first appearance is not in the context of a physical context at all, 
but in the study of the geometry of surfaces with constant negative Gaussian curvature. Physical 
problems for which it has been used include: Josephson junction transmission lines, dislocation in 
crystals, propagation in ferromagnetic materials of waves carrying rotations in the magnetization 
direction, etc.^° Mathematically, it is a very interesting because it is one of the few physically 
important nonlinear partial diff'erential equations that can be solved explicitly (by a 
technique known as Inverse Scattering, which we will not describe here). 

An important consequence of equation (3.12) exact solvability, is that it possesses particle- like 
solutions, known as kinks, anti-kinks, and breathers. These are localized traveling distur- 
bances, which preserve their identity when they interact. In fact, the only effect of an interaction 
is a phase shift in the particle positions after the interaction: effectively, the "particles" approach 
each other, stay together briefly while they interact (this causes the "phase shift") and then depart, 
preserving their identities and original velocities. This can all be shown analytically, but here we 
will only illustrate the process, using some computational examples. 



^^Fbr reviews see: 

A. C. Scott, 1970, Active and Nonlinear Wave Propagation in Electronics, Wiley Interscience, New York (page 250). 
Barone, A. F. Esposito, C. J. Magee, and A. C. Scott, 1971, Theory and Applications of the Sine Gordon Equation, 
Rivista del Nuovo Cimento vol. 1, pp. 227-267. 
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The first step is to present analytical expressions for the various particle-like solutions of 
equation (3.12). These turn out to be relatively simple to write. 

Example 3.1 Kinks and Anti-Kinks. 

Equation (3.12) has some interesting solutions, that correspond to giving the pendulums a full 27r 
twist (e.g.: take one end pendulum, and give it a full 'In rotation). This generates a 2ti twist wave 
that propagates along the pendulum chain. These waves are known as kinks or anti-kinks (depending 
on the sign of the rotation), and can be written explicitly. In fact, they are steady wave solutions,^^ 
for which the equation reduces to an O.D.E., which can be explicitly solved. 



be a moving 



Let I — 1 < c < 1 1 be a constant (kink, or anti-kink speed), and let z = {x — ct — Xq) 

coordinate, where the solution is steady — the "twist" will be centered at x = ct + Xq, where Xq is 
the position at time t = 0. Then the kink solution is given by 



6 = 2 arccos ( 2z/p _^ i j ~ ^ arctan ( exp (""^11 ' (3.13) 



where 



c 



is the kink width. This solution represents a propagating clock-wise 27r rotation, 
from 9 = 2rmr as x ^ — oo (where m is an integer) to 9 = 2 (m — 1)% as x ^ oo, with most of the 
rotation concentrated in a region of width 0{f3) near x = ct + xq. The parameter c is determined 
(for example) by how fast the initial twist is introduced when the kink is generated. 

We note now that: 

• From (3.13) it follows that 9t = —c9x = sin ■ Using this, it is easy to show that (3.13) 
is a solution of equation (3.12). 

• The Sine-Gordon equation is the simplest of a "class" of models proposed for nuclear inter- 
actions. In this interpretation, the kinks are nuclear particles. Since (in the non-dimensional 
version (3.12)) the speed of light is 1, the restriction —l<c<lis the relativistic restriction, 
and the factor /5 incorporates the usual relativistic contraction. 

The anti-kink solution follows by replacing x — ^ —x and t ^ —t in (3.13). It corresponds to a 
propagating counter- clock-wise 27r rotation, and it is given by 

/l_e22//3\ / /^\\ 

9 = 2 arccos r— = 4 arctan exp — • (3.14) 

The kinks and anti-kinks are very non-linear solutions. Thus, it is of some interest to study how 
they interact with each other. Because they are very localized solutions (non-trivial only in a small 
region), when their centers are far enough they can be added. Thus, numerically it is rather easy to 
study their interactions, by setting up initial conditions that correspond to kinks and anti-kinks far 
enough that they do not initially interact. Then they are followed until they collide. In the lectures 
the results of numerical experiments of this type will be shown (the numerical method used in the 
experiments is is a "pseudo-spectral" method). 



^^Solutions of the form 6 = 6{x — ct), where c is a constant: the speed of propagation. 
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Example 3.2 Breathers. 

A different kind of interesting solution is provided by the "breathers" — which we handle next. A 
breather is a wave-package kind of solution (an oscillatory wave, with an envelope that limits 
the wave to reside in a bounded region of space. These solutions vanish (exponentially) as x ^ ±cx). 
This last property allows for easy numerical simulations of interactions of breathers (and kinks). 
One can setup initial conditions corresponding to the interaction of as many kinks and/or breathers 
as one may wish (limited only be the numerical resolution of the computation), simply by separating 
them in space. 

A breather solution is characterized by two arbitrary constants —l<d,V<l. Then define 



A 
B 
C 

P 



d/VT^ 

1/VT 



v. 



d^. 



CB{Vx-t + to) 
dB {x — V t — Xq) . 



(3.15) 



Q = y4 sin (p)/ cosh (g) , 

where xq and to are constants, centering the envelope and the phase, respectively. Notice that the 
partial derivatives of Q (with respect to p and q) are given by 



Qp = A cos(p)/ cosh(g) and Qg = —Q tanh(g) 
The breather solution (and its time derivative) is then given by: 

9=4 arctan((3) , 

9t = -A{l + Q^){CBQp + dBVQg). 



(3.16) 



(3.17) 



The breather solution is a wave-package type of solution, with the phase controlled by p, and the 
envelope (causing the exponential vanishing of the solution) by q). The wave-package details are 
given by: 



speed 


. . . . Cp 




period 


T 


= 2 7r/(SC), 


wave-length 


\ 


= 2n/{BCV), 


speed 


Ce 


= V, 


width 


Ae 


= 27r/{dB), 



Phase. 



Envelope. 



(3.18) 



(3.19) 



Notice that, while the phase moves faster than the speed of "light" (i.e.: 1), the envelope always 
moves with a speed —1 < V < 1, a,nd has width proportional to a/I — V^. 

Finally, in case you are familiar with the notion o/ group speed, notice that (for the linearized Sine- 
Gordon equation: 9u — 9xx -\- 9 = Q) we have: (group speed) = l/(phase speed) — which is exactly 
the relationship satisfied by Cg = V and Cp = 1/V for a breather. This is because, for \x\ large, the 
breathers must satisfy the linearized equation. Thus the envelope must move at the group velocity 
corresponding to the oscillations wave-length. 
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Remark 3.5 Pseudo-spectral Numerical Method for the Sine-Gordon Equation. 

Here we will give a rough idea of a numerical method that can be used to solve the Sine-Gordon 
equation. This remark will only make sense to you if you have some familiarity with Fourier Series 
for periodic functions. 

The basic idea in spectral methods is that the numerical differentiation of a (smooth) periodic func- 
tions can be done much more efficiently (and accurately) on the "Fourier Side" — since there it 
amounts to term by term multiplication of the n*'* Fourier coefficient by in. On the other hand, 
non-linear operations (such as calculating the square, point by point, of the solution) can be done 
efficiently on the "Physical Side". 

Thus, in a numerical computation using a pseudo-spectral method, all the operations involving taking 
derivatives are done using the Fourier Side, while all the non-linear operations are done directly 
on the numerical solution. The back-and-forth calculation of Fourier Series and their inverses is 
carried by the FFT (Fast Fourier Transform) algorithm — which is a very efficient algorithm for 
doing Fourier calculations. 

Unfortunately, a naive implementation of a spectral scheme to solve the Sine-Gordon equation would 
require periodic in space, solutions. But we need to be able to solve for solutions that are mod-27r 
periodic (such as the kinks and anti-kinks), since the solutions to the equation are angles. Thus, 
we need to get around this problem. 

In a naive implementation of a spectral method, we would write the equation as 



Ut = v, 

Vt = — sinu . 



(3.20) 



where u = 9 and v = 9t. Next we would discretize space using a periodic uniform mesh (with a large 
enough period), and would evaluate the right hand side using FFT's to calculate derivatives. This 
would reduce the P.D.E. to some large O.D.E., involving all the values of the solution (and its time 
derivative) at the nodes in the space grid. This O.D.E. could then be solved using a standard O.D.E. 
solver — say, ode45 in MatLab. 

In order to use the idea above in a way that allows us to solve the equation with mod-'lir periodicity 
in space, we need to be able to evaluate the derivative u^x in a way that ignores jumps by multiples 
of 21: in u. The following trick works in doing this: 



Introduce 



U = e*". 



Then 



(Ux?-UUxx 

Uxx = i ^ ' ^2 (3-21) 

gives a formula for Uxx that ignores 27r jumps in u. Warning: In the actual implementation one 
must use 

[ {Uxf-UUxx \ 

Uxx = -imag I — I 

to avoid small imaginary parts in the answer (caused by numerical errors). 



4 Suggested problems. 



A list of suggested problems that go along with these notes follow: 
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1. Check the derivation of the system of equations (2.20). 

2. Derive the continuum equation in (2.21). 

3. Look at the end of section 2, under the title "General Motion: String and Rods". Derive 
continuum equations describing the motion (in the plane) of a string without constraints. 

4. Look at the end of section 2, under the title "General Motion: String and Rods". Add bending 
springs to the model, and derive continuum equations describing the motion (in the plane) of 
a rod without constraints. 

5. Do the check stated in remark 3.2. 

6. Answer the question in remark 3.3. 

7. Do the dimensions check stated below equation (3.11). 

8. Answer the question in remark 3.4. 

9. Show that (3.13) is a solution (there is a hint about how to do this a few lines below the 
equation). 

10. Use a computer to plot the solution in (3.13), as a function of z, for a few choices of c. 

11. Show that (3.17) is a solution. 

12. Use a computer to plot the solution in (3.17), as a function oi x, for various times and choices 
of parameters. 

13. Implement a numerical code to calculate interactions of kinks, breathers, etc., using the ideas 
sketched in remark 3.5. 



